Theoretical analysis of super— Bloch oscillations 



K. Kudo 1 and T.S. Monteiro 2 

1 Division of Advanced Sciences, Ochadai Academic Production, 
Ochanomizu University, 2-1-1 Ohtsuka, Bunkyo-ku, Tokyo 112-8610, Japan 
2 Department of Physics and Astronomy, University College London, 
Gower Street, London WC1E 6BT, United Kingdom 
(Dated: January 11, 2013) 

Several recent studies have investigated the dynamics of cold atoms in optical lattices subject 
to AC forcing; the theoretically predicted renormalization of the tunneling amplitudes has been 
verified experimentally. Recent observations include global motion of the atom cloud, such as giant 
"Super-Bloch Oscillations" (SBOs). We show that, in order to understand unexplained features of 
SBOs, in addition to the renormalization of the tunneling, a new and important phase correction 
must be included. For Fermionic systems with strong attractive interactions, one may engineer 
different types of collisions and recollisions between bound-pairs and unpaired atoms. 
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I. INTRODUCTION 

Recent experiments on cold atoms in optical lattices 
subject to time-periodic perturbations [lHy] have pro- 
vided clean realizations of a range of different types of 
AC control of coherent matter waves proposed by earlier 
theoretical studies 043 ■ Even neglecting the effects of 
inter-particle interactions, these studies identified non- 
trivial dynamical effects. In particular, for an oscillating 
potential of strength F and frequency uj, it was found 
that the tunneling amplitudes J of the driven atoms take 
an effective, renormalized, value: 
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where d is the lattice constant and J n denotes an ordi- 
nary Bessel function and n = 0, 1,2... Values of n > 
arise if an additional static linear (in x, the position) 
field Fox is applied, satisfying a resonance condition 
F d = nhjj. The inter-site transport is completely sup- 
pressed at parameters corresponding to the zeros of the 
Bessel functions. 

The oscillating potential can be implemented with ul- 
tracold atoms in shaken optical lattices, which yield po- 
tentials of the form Hp(t) = Fxsin(ujt + 0), where 
4> is a phase which can be controlled. Thus, Eq.([T]) 
was demonstrated and investigated experimentally in 
Refs. [l|, The observed suppression of transport is 
variously termed Dynamic Localization (DL) or Coher- 
ent Destruction of Tunneling (CDT), depending on the 
transport regime 1111 • They are both effects of the 
one-particle dynamics; comparatively less work has been 
done in the strongly interacting regime though, for ex- 
ample, in 0, Il3l - |l5j | they have investigated interactions- 
driven effects. Recent studies consider even triangular, 
shaken lattices [l6j . 

The experiments in Refs. (HQ investigated the spread- 
ing or local tunneling of an atomic wavepacket, without 
global motion. But an earlier theoretical study had 



proposed also the possibility of global transport of the 
atomic wavepacket in the presence of the additional static 
field and assuming i*o^ — ^ li-Sy n = 1). Recent ex- 
periments using a static field [J, Q were able to realize 
both directed motion as well as large oscillations, occur- 
ring over hundreds of sites, which were termed "Super- 
Bloch Oscillations" (SBOs) %. These SBOs were ana- 
lyzed in Ref. [l7| . including the effects of weak interac- 
tions (mean-field regime). The general assumption was 
that the group velocity in these cases v g oc ^Ti(F u ), i.e., 
the dependence of the dynamics on the oscillating field 
F^ = Fd/(Hu)), is entirely contained in the Bessel func- 
tion argument. 

We show here that, in order to explain the new exper- 
iments, a phase correction F^coscj) — n{<j) + £) must — 
additionally — be considered cither in an effective disper- 
sion relation or in the average group velocities. We show 
that inclusion of the appropriate terms can account for 
several unexplained experimental features in the exper- 
iments [3, S|- A notable example is a field dependent 
shift observed in the phase of SBOs, not accounted for 
by the usual analysis. Other features include the differ- 
ent field dependence of the speed of directed motion, on 
resonance, relative to the amplitude of SBOs. 

There is also currently much interest in pairing phe- 
nomena, motivated by many ground-breaking experi- 
ments with ultracold Fermionic atoms in optical lattices 
[H,[lj|- We have previously shown [2(| that in the case of 
attractive interactions, where the atoms can form bound- 
pairs, the second-order tunneling mechanism implies a 
different renormalization by the field and thus develop a 
global velocity relative to unpaired atoms. We investi- 
gate here for the first time the dynamics of pairs in the 
SBOs regime. We find the paired/unpaired components 
can be made to re-collide repeatedly; the novelty here is 
that the new f^-dependent phase term enables control 
of re-collisions: the two components can either touch and 
reflect or can be forced through each other. Experimen- 
tal demonstration is a matter of combining the techniques 
used to study of dynamics of bound pairs fl9| with the 
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AC driving in Refs. [iHa]; the re-collisions in two com- 
ponent ultracold gases suggest the possibility of other 
applications in cold chemistry. 

The central objective of the present work is the effect of 
the phase of the AC driving acting on an initially undis- 
turbed atomic wavepacket. Thus, the phase <j> must be 
well defined over every cycle, and the driving field must 
be switched on initially over a time that is small relative 
to the driving period. But provided this constraint is 
satisfied, the switch-on protocol is immaterial. The ob- 
served effects are never due to a discontinuous "switch- 
on" or an abrupt jump in amplitude of the potential at 
t = 0; in fact the strongest experimental effects occur for 
<f> = 0, 7r where the potential ramps up smoothly from 
zero at t = (i.e., a pure sinwi drive). The authors 
of Ref. studied the effect of adiabatically switching 
on AC driving over many driving periods on the Floquet 
quasi-energy bands, where the phase is thus of no conse- 
quence. The authors of Ref. [22J considered slow linear 
ramping of the driving f(t) = Ktsin(ujt + cj>) by means 
of a perturbative treatment valid for high lo. For a ramp 
period lasting an integer number of periods t = nT, a net 
ratchet current was obtained. Other forms of driving like 
amplitude modulated lattices that do not show dynamic 
localization can also produce directed motion [23[. 

In Section II below, we review key aspects of the the- 
ory of renormalization of tunneling for atoms subject to 
AC forcing. In Section III we show that the phase corre- 
sponds to directed motion and compare with other sys- 
tems. In Section IV we analyze Super-Bloch Oscillations 
and obtain an expression which we show explains impor- 
tant features of recent experiments. 



II. RENORMALIZATION OF TUNNELING 

As in Refs. [5. [Tl|. we consider the dynamics in a spa- 
tially periodic potential, subject to an additional time- 
periodic driving term. The total Hamiltonian is H(t) = 
Ho + Hp(t). Here, H corresponds to the non-interacting 
limit of a variety of Hamiltonians with nearest-neighbor 
hopping (Hubbard, Bose-Hubbard, magnons in Heisen- 
berg spin chains, etc.). It represents any spatially pe- 
riodic potential characterized by energy eigenfunctions 
fymk-, with band index m and wavenumber k, thus 
Ho4> m k = E m {k)4> m }~. We restrict our one-particle prob- 
lem, i.e. 



so this corresponds to the group velocity 
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where cj and Cj are the creation and annihilation oper- 
ators of a Fermion or Boson, to the lowest band m = 1; 
taking E(k) = E m —i(k), the energy dispersion 
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Assuming that the external driving is linear in position, 
we have H F (t) = —f(t)x where f(t) = F a +F sm(ujt + (f)); 
in general it comprises both a static field and a sinu- 
soidally oscillating field with an arbitrary phase <f>. The 
result of the driving is a time-dependent wavenumber: 



qk{t) = k + JiJo dTf(r) ' 



(5) 



The stationary states of the system are its Floquet states, 
the analogues of Bloch waves in a temporally periodic 
system. They are given by 



ip k (x,t) = u k (x,t)exp 
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where Uk(x, t) — Uk(x, t 4- T), and the period T = 2ir/ui. 
The non-periodic phase term is characterized by the 
quasienergy e(k). The evolution of a wavepacket pro- 
jected onto its Floquet states is fully determined: the 
quasienergies play a role entirely analogous to the energy 
eigenvalues of a time- independent system; the (period- 
averaged) group velocity v g of a wavepacket is given from 
their dispersion: 
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evaluated at the appropriate initial momentum, ko, in 
analogy to Eq. (j4|) for the undriven system. Below we 
take h — 1, which implies F u = Fd/ui. 

The presence of the static linear field, in general, de- 
stroys the band dynamics; however, here we consider 
the so-called resonant driving case, for which Fad = nuj, 
(where the driving compensates for the energy offset be- 
tween neighboring wells in the lattice, restoring the band 
structure). In this case, it can be shown that an ef- 
fective quasienergy dispersion is obtained from the energy 
dispersion by a period-averaging, over one oscillation: 



e(k) = ^J E(q k (t))dt. 



(8) 



We begin by considering the case n = : in the first 
experimental studies on Dynamic Localization, Fq = 
and thus the static field was absent. In previous theoret- 
ical studies (Ref. Q) a driving term of form sin(wi + (f>) 
was considered, but because of the particular objectives 
of that work, the effects of the phase <j> were not re- 
tained. In that case, evaluating the integral in Eq. (jHJl, 
the well-known renormalization expression (see also de- 
tailed derivation in Ref. [31) was obtained: 



E(k) = — J cos kd, 



(3) 



e(k) = - J J (F U ) cos kd. 



(9) 
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It can be seen that the tunneling amplitude is multi- 
plied by a Bessel function. Now the hopping can be com- 
pletely suppressed at the zeros of the zero-th order Bessel 
function. This process was first demonstrated in Ref. [l| . 
However, for the typical initial wavefunction correspond- 
ing to a zero-momentum ultracold atom cloud, for exam- 
ple, a Gaussian sharply peaked about k — ko = 0, Eq. ([7]) 
indicates an average group- velocity v g — 0. 



III. DIRECTED MOTION 

The situation is quite different if the inte gra l Eq. ([8]) is 
evaluated without disregarding cf>. In Ref. [20j we found 
that even for the case <p = 0, a momentum shift in the 
effective dispersion results. For the case of general </>, the 
effective dispersion relation becomes: 



e(k) = — JJo(F u ) cos [kd + cos i 
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where one sees that there is an F u cos <f> shift, representing 
the average momentum over one cycle. This is no longer 
equal to zero. The result of this field-dependent shift is to 
introduce directed motion, at constant (period-averaged) 
group velocity [20(. 

Atoms in shaken lattices experience a homogeneous 
(position independent) force. Independently, theoretical 
[24j and experimental [25| studies considered the role of 
phase jumps in driven traps. These have inhomogeneous 
forces; for a driving field, V(x, t) = V'(x) sin(o;i+</>) clas- 
sically, the phase effects a position-dependent momentum 
shift AP(x) oc V'(x). The position dependence of V'(x) 
strongly couples the phase to the dynamics. AP(x) can 
generate larger /smaller amplitude oscillations in the trap 
and has been proposed as a means to control the kinetic 
energy of the atoms oscillating in the trap. 

Despite the different dynamics, for all the above sys- 
tems, the effect of the phase vanishes for <f> — ir/2. 

In Eq. (|10p , the F u cos <fi phase shift coincides with the 
value of the lower bound of the integrand in Eq. ([5]) (here 
evaluated at t = 0). This might lead one to conclude that 
the shift arises from the abrupt jump in the amplitude 
of the driving potential at t ~ 0. But this would be 
a misapprehension; the key physical significance of the 
shift is that it is the average momentum over each cycle, 
since (p(t))T = F u cos</> takes the same form. Below we 
find that the strongest experimental effects occur for (f> — 
0, 7r where the driving grows linearly from zero for t ~ 
(a pure sinwi drive). The phase must be well defined 
over every cycle, so the switch-on time to should satisfy 
the condition to -C T. But provided this constraint is 
satisfied, the switch-on protocol is immaterial. 

We now analyze recently discovered large scale oscilla- 
tions for cold atoms in optical lattices. Although they are 
phase-dependent, we note that they are quite different 
from the atoms in traps: the dynamics are independent 
of initial position (the wavepacket is delocalized over sev- 
eral wells). We show also that the SBOs rely crucially 
on a specific quantum resonance (more precisely a slight 



detuning from it). We show below that their amplitude 
does not depend on 4> (unlike oscillations in traps); only 
their phase does. 



IV. SUPER-BLOCH OSCILLATIONS 

Below we also consider the case of non-zero Fq as well 
as slight detuning for which F^d = (n + S)uj, with 5 <C 1, 
associated with SBOs, for which the above relation still 
holds. 

In order to calculate (period-averaged) group veloci- 
ties, we first evaluate Eq. ([5]): 

q k (t) = k +[{n + 5)u)t - F u cos(ujt + <j)) + F^ cos <j)] /d, 

(11) 

(assuming qk(t = 0) = k) then substitute the result in 
Eq. ©. However, for the slight-detuning case 5^0, 
we assume that the time-dependence due to the Suit/d 
remains negligible over one period T. Thus, we take it 
out of the integral and Eq. §8$ becomes (see Appendix [X] 
for the detailed derivation): 

e(k) ~ —JJn{Fu,) cos kd + 5ujt + F^ cos0 — n(<p + — ) 

(12) 

The above represents an effective dispersion relation, but 
which oscillates slowly in time with a period Tsbo — 
2tt/(Slu) > T B , where T B oc 1/F is the Bloch period. 
They correspond to the SBOs investigated by Refs. 0, H, 
d, OS]- Even at resonance 5uj = 0, Eq. (fT2|) differs from 
previous expressions by the phase-shifts F u cos (f> — n((j> + 

Evaluating Eq. (TT^|) for the n — 1 case, we obtain 
de 

v g = — ~ - Jdcos (kd + F u cos <j> — <f> + Scot) . 

(13) 

However, experiments [5( measure the center-of-mass po- 
sition x(t) = J Q v g (t')dt', where the inter-site spacing 
d ~ 0.533 /im is included if x(t) is obtained in /im. In 
that case we obtain 

x{t) ~ —^-Ji{F w ) [sm(K F + 8ut) - smK F ] , (14) 

OUJ 

where Kp = kd + F^ cos (f) — <fi = F u cos (j> — 4>, since in 
experimental situations of interest here k ~ initially. 

Equation. (fl"4"|) can largely account for the complex de- 
pendence of experimental results on <f>. For instance, it 
was noted in Ref. Q that the experimental phase of 
the SBOs depends on the sign of S, for k = 0; this 
would not be expected without the F^ cos <f> — 4> shifts, 
since cos(5ujt) = cos(-Sujt) [see Eq. (|12))]: it was also 
noted that the SBO amplitudes scale as l/(<5w); also, 
for <j) = 0,7r, the SBO amplitude was close to a maxi- 
mum for F u = 1.52 ~ 7r/2. In contrast, we see that for 
6 = (directed motion) the peak velocity occurs wher- 
ever cos(F UJ )Ji(F LU ) is a maximum (i.e., at F u ~ 1). The 
directed motion is in fact almost zero for F^ = 1.52, the 
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point where the SBOs were near their maximum, but 
the highest experimental directed motion was given for 
F u ^ 1. 

However, the most interesting experimental feature 
predicted by our Eq. (IT4l) is that the SBOs begin with 
a field-dependent phase, a feature surprisingly evident 
even in measurements not looking for this behavior. Fig- 
ure [JJ demonstrates this for F u = 0.15, 1.52, tt/2 and tt, 
in the non- interacting Hubbard model. Figures [TJ (c) and 
(d) show that Eq. ([TT]) reproduces quite well the experi- 
mental values of Ref. @ especially for small times. The 
disagreement with experimental data in large times can 
be attributed mostly to interactions. The graph shows 
clearly the displacement of the first maximum, seen in 
the experiment as well as the order of magnitude vari- 
ation in amplitude. Such a field- dependent shift is also 
apparent in the results of Ref. Figures [TJ (e) and 
(f) contrast = tt/2 with F^ = tt. For the former, the 
atomic wavepacket begins with near-zero speed and grad- 
ually accelerates, while for tt it leaves with maximal speed 
(i.e., it acquires the speed on a timescale Tb <C Tsbo) and 
returns to the original position with maximal speed. As 
we see below, this has interesting physical implications 
if, when the wavepacket returns to its original position, 
it re-collides with an immobile component. 



V. BOUND-PAIRS 

We now consider the effect of the SBO phases in a 
regime of attractive interactions and Fermionic atoms, 
such as investigated in Refs. [TH, In that case, we 
model the dynamics by the many-body Hubbard Hamil- 
tonian: 
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where cj a and c Jj(T are the creation and annihilation op- 
erators of a Fermion and rij. a = c^^Cj^. This supports 
Bound-Pair (BP) states, lower in energy by \U\ relative 
to the unpaired states; as for the unpaired atoms, we 
consider their motion in the lowest band only, with dis- 
persion relation p6l |: 

J2 J2 

E(k 1 ,k 2 ) = -— cospi +k 2 )d] = -— cos2/id, (16) 

where k is the center-of-mass momentum. Similarly, in 
the many-body Hamiltonian, the driving term: 



Hi{t) = - [F Q +Fsm(ujt 



j 



(17) 



where j denotes the site index and hj is the number oc- 
cupancy ( = 2 for a BP; and = 1 for a site j occupied 
by either an unpaired spin- up or spin-down atom) . Thus 
Hf p {t) = 2£/ 1 unpaircd (i), i.e., the magnitude of the driv- 
ing term is doubled for the BP. 
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FIG. 1: (Color online) (a) and (b) One-particle solutions 
of the Hubbard Hamiltonian for L — 30 lattice sites, show- 
ing the field-dependence of the phase of SBOs. Upper/lower 
panel correspond to F u = 1.52 and 0.15; both have <f> = tt. For 
illustration purposes, the theoretical amplitudes were equal- 
ized by choosing J- values which equalize J e g ~ JJi(F u ). (c) 
and (b) Shows that good agreement is obtained between 
experiment and the Hubbard numerics, as well as the ana- 
lytical formula: Eq. (I14|l (dashed lines) reproduces well both 
the amplitude and phase of SBO experimental data (symbols) 
of Ref. Q, using 8oj = 2tt/1000 (s-IHz detuning) 4> = tt 
and Jd/(h5uj) = 90 fim. This implies J/h « 1.06 ms . 
(e) and (f) show one-particle Hubbard numerics contrasting 
Fu = tt/2 (upper panel) with n (lower panel), parameters 
used below for collisions with bound-pair states. For F u = tt, 
the wavepacket leaves at high speed and will return (and rec- 
ollide with any component) at the original position with high 
speed; for = n/2 the wavepacket starts with zero speed 
and returns with zero speed. 



Repeating the arguments leading to Eq. (|12|) . we can 
show that for the regime of the experiments (n = 1 and 
(f> = 0) 0, Hi , the BP group velocity is 



..bp 



J 2 d 

-—sin (2nd +2F u + 26u>t) J 2 (2F UJ ). (18) 



Thus, the order of the Bessel, its argument, and the ad- 
ditional phase are all doubled; so the zeros and maxima 
of the BP motion occur at different fields from the un- 
paired atoms. However, in the presence of the driving, 
neither the BPs nor unpaired states remain eigenstates 
of the system. For |J7| > 2(Fo + F), though, the energy 
gap between BPs and unpaired atoms cannot be closed 
by the fields, and there is relatively little mixing. We 
tested this assumption for simulations for lattices with 
L ~ 30-50 sites and N = 2,3, and 4 atoms. Numer- 
ics are presently unfortunately limited by N = 4 and 
L ~ 30; since N = 4 includes the residual BP-BP inter- 
actions (and BP's interact like hard-core bosons), it is 
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FIG. 2: (Color online) Solution of the Hubbard Hamiltonian 
dynamics of a mixture of Bound-Pairs and unpaired atoms in 
the SBOs regime for L — 30 or 50 lattice sites. The dynamics 
are qualitatively similar for N = 2-4 atoms; these calculations 
use N = 3, with two spin-up and one spin-down atoms. Since 
J 2 d/U <C 1, the bound-pair component is almost static and 
appears as a bright horizontal line; the free atoms (l/3rd of 
the total probability, initially) perform SBOs around them, 
(a) and (b) correspond to U = —10, F u = 7r. The unpaired 
atoms separate then recollide with BPs at maximum speed, 
then attempt to cross the static BP component. At low energy 
and large effective h oc 1/L ((a) J = 0.6, L = 30), there is 
a mixture of transmission and reflection; At higher energy 
and small fi ((b) J = 1, L = 50), complete reflection and a 
series half SBOs predominate, (c) and (d) are for U = — 7, 
F u — tt/2. For F w — tt/2, the wavepacket starts at zero speed 
and returns with zero speed for both L = 30, J = 0.5 ((c)) 
and L = 50, J = 1 ((d)); there is relatively little transmission. 



attempt to "push-through" the BPs, leading in general to 
a degree of beam-splitting, which depends on the kinetic 
energy (i.e., J) and effective h (oc 1/L). The F u = tt/2 
atoms, on the other hand, return with zero velocity and 
simply turn around. 

Our analysis on BP dynamics is valid only for \U\ 3> 
lo ~ J . In this regime, bound pairs remain in the BP 
band, and thus Eqs. (fT6|) and (|T8| are likely to remain 
valid. If u> 3> |Z7|, new types of tunneling might be found. 
Such a high-frequency regime warrants further investiga- 
tion. 



VI. CONCLUSION 

We have shown a phase-correction in an effective dis- 
persion relation is essential for full understanding of cur- 
rent experiments on transport with AC forcing which 
showed Super-Bloch Oscillations. 

In addition, we investigate Super-Bloch Oscillations 
for systems with bound pairs. We show that we can con- 
trol collisions and re-collisions between unpaired atoms 
and bound pairs. This has potential implications for 
studies of AC control of two-component condensates, in- 
cluding Fermionic systems and molecular condensates, of 
relevance in cold chemistry. 

We acknowledge helpful comments or discussions 
with A. Eckardt, M. Holthaus, C. Weiss C. CrefEeld, 
A. Kolovsky, G. Ferrari, E. Arimondo and O. Morsch. 
We are extremely grateful to E. Haller and C. Nagerl 
for the experimental data in Fig. [TJ This work is partly 
supported by KAKENHI(21740289). 



expected that the essential physics is included. 

Nevertheless, we take a conservative approach and re- 
port below only on dynamics that are qualitatively insen- 
sitive to whether N = 2 or N — Thus, we consider 
the large U regime where BP and unpaired dynamics 
is separated and break-up rates of the BPs are modest: 
for \U\ < 2(F + F), the pairing rate falls rapidly and 
can oscillate in time. A further advantage of large \U\ is 
that the BPs are essentially static, thus reducing another 
source of uncertainty. 

In Fig. [2] we show the results of calculations slightly 
above the threshold where the field can strongly cou- 
ple the Bound-Pair states and free atoms, and there is 
already a certain degree of interaction in evidence. It 
is seen that the BP component is essentially immobile, 
while the unpaired atoms perform large-scale oscillations; 
the unpaired atoms periodically return and re-collide 
with the BP packet. However, the F u = n trajectories 



Appendix A: Effective dispersion relation 

In order to derive Eq. ([T2"]) from Eq. ([TTj) . we assume 
that the Super-Bloch period represents a completely 
different timescale from the much faster Bloch period. 
Thus, we introduce two timescales, a rapid time t for the 
time-average over a period and a slow time t' = Suit which 
describes the slow dynamics. Equation (|12p is rewritten 
as 

q k (t) = k+{t' + F,, cos fi + nujt-Fu cos(ujt + ((>)}/ d. (Al) 

Since 6 <C 1, so t' <C t and the change due to t' during 
one period T is negligible. Substituting (|A1[) into ©, we 
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have where r = tot + 4>. Evaluating both integrals above, one 

obtains Eq. (IT21) : 

J f 

e(k) = — — / dt cos [kd + t' + cos <j) 
1 Jo 

+nojt — cos(ujt + 4>)] 

= -— [ drcos[kd + t' + F U cose/)- ncj) ^ ~ - J ^ F ") 

27r J<p 

+71T — F u cos t] 
~ — cos(kd + t' + F w cos — n0) 

Z7T 



X COS 



kd + 6ut + F u cos <(>-n((f)+^) (A3) 



x 



/ dr cos [nr — F u cos r] 
■/ o 



J f The Bessel function in Eq. (TTzT) is obtained from the 

2^: sm (fc^ + * + F w cos — n(/)) (non-zero) integrals; the additional term mr/2 in the 

p2ir phase arises because the first integral vanishes for n odd, 

x / drsin[nr — F u cost], (A2) while the second integral vanishes for n even. 
Jo 
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